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Level density and radiative strength functions are obtained from the analysis of two- 
step cascades intensities following the thermal neutrons capture. The data on level density 
are approximated by the sum of the partial level densities corresponding to n quasi-particles 
excitation. The most probable values of the collective enhancement factor of the level density 
are found together with the thresholds of the next Cooper nucleons pair breaking. These 
data allow one to calculate the level density of practically any nucleus in given spin window 
in the framework of model concepts, taking into account all known nuclear excitation types. 
The presence of an approximation results discrepancy with theoretical statements specifies 
the necessity of rather essentially developing the level density models. It also indicates the 
possibilities to obtain the essentially new information on nucleon correlation functions of the 
excited nucleus from the experiment. 

1 Introduction 

First of all, obtaining data on the level density p and radiative strength functions k emission 
of nuclear reaction products is a way to universally test nuclear models. Their accuracy is 
quite unequivocally limited by the experimental data quality used for both creation of new 
modeling representation, and for its parameterization. Naturally, the experimental data used 
for this purpose should be the most reliable. At present, for evaluating of [1] neutron cross 
sections it is required that there is a high reliability of modeling representations on the level 
density in a wide nucleus excitation energy interval. It is also important for describing [2] 
their fission process. That means they are needed for quite practical applications of nuclear 
physics. Thus, materials' and nuclear power cycle fuel's profitability and safety directly 
depend on the accuracy of their cross sections evaluation. 

The history of developing the techniques for experimental determination of p covers 
several decades, but the significant part of the experimental data obtained up to now has 
too large systematical error, which is essentially unavoidable within the framework of existing 
experimental analysis algorithms. It is connected to extraction peculiarity of the level density 
and radiative strength functions from experimentally measured spectra. This regards the 
situation, when the insufficient resolution of nuclear reaction's products, provided by the 
existing spectrometers, does not allow one to resolve the individual levels of a nucleus. In 
wide energy intervals, where the nuclear reactions' products are registered, the intensity of 
spectra simultaneously depends on both p and k. Errors transfer results in a large increase 
of systematical errors in the found values. It always occurs when there are systematical 
experimental errors of the measured spectra and the strong correlation of parameters. 
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Currently, there is a set of various advanced modeling representations [3] for predicting 
both the level density and the radiative strength functions of the primary gamma-transitions 
following compound states decay. However, the accuracy of k values, calculated with their 
help, is unknown for the low-energy primary gamma-transitions following neutron resonance 
decay Therefore, contemporary techniques, in which parameters of a nucleus are extracted 
from the measured spectra of the nuclear reactions gamma-rays, can provide no high accuracy 
for either level density or the radiative strength functions. It is impossible if the modeling 
representations on one of these values are used to determine the other. 

These circumstances demand an obligatory revealing of all of the most essential sys- 
tematical errors of p and k certain values by comparing any model representation with the 
experiment. 

2 Systematical errors of the level density and radia- 
tive strength functions in the experiments of various 
types 

By present time four different techniques to determine functional dependence of the level 
density are realized. They are based on the following experimental data: 

1. Nuclear evaporation spectra in reactions of various types; 

2. Full spectra of gamma-transitions for various energies of the decaying levels; 

3. Intensities of the two-step cascades appearing at the thermal neutrons capture; 

4. Distribution of the random intensity values of individual energetically resolved two- 
step cascades with the energy of their intermediate level Ei < 0.5B n . 

Here, it is necessary to add a completely independent technique of [4] to determine the 
deviation sign of the experimental level density from any model prediction. The corre- 
sponding analysis uses the data on random fluctuations of the intensity ratios of secondary 
gamma-transitions at the resonance neutrons capture. Its results obviously contradict the 
conclusions of the 1-st and 2-nd techniques and qualitatively confirm the basic of the 3-rd 
and 4-th methods that the level density is smaller than that given by the present standard 
models. 

As it is stated above, all specified techniques always result in p and k values having 
unknown systematical errors. The most exact of their estimation can be obtained only by 
comparing their values, determined in various and complementary independent experiments. 
The comparison quite uniquely shows that the level density in the region of the half of neutron 
binding energy, found within the framework of the first two techniques is ~ 3 — 5 times (in 
a number of nuclei - and more) larger than that determined according to the 3-rd and 4- 
th algorithms. At lower and higher excitation energies this discrepancy decreases because 
densities of resolved low-lying levels and of neutron resonances are always normalized to the 
same values in techniques 1-3. 

Thus, during development and parameterization of p models, results of application of the 
level density, obtained by using of the listed above techniques, may bring to rather fallacious 
ideas about this parameter of a nucleus. 

1. Practically, all theoretical models of level density [3] are completely developed on the 
base of the experimental data on p, obtained similarly to [5] from evaporation spectra of 
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nucleons or light nuclei in various nuclear reactions. In order to determine the experimental 
level density, the model-set values of transmission coefficient T of a nucleus surface have 
been used for the emitted nucleons or light nuclei. Up to the present the analytical method 
of an optical model of a nucleus is used for their determination, providing that there is no 
determinate choice of nuclear potential used for calculation. It is necessary to take into ac- 
count, that 283 various parameterizations of optical potential were developed to describe the 
neutron-nucleus interaction process; 101 - for proton, correspondingly [3]. Such a quantity 
of entrance factors needed for calculation of the transmission coefficient in reaction (n,p) [5] 
serves as a serious reason for assuming that there are very big errors in extraction of level 
density by means of the first technique. Naturally, determination of level density with the 
minimal error requires one to calculate T values within the framework of nuclear models 
with the guaranteed accuracy which exceeds that needed for p determination. The latter, 
by magnitude scale, may be equaled to the modern accuracy of observation of the spectrum 
intensity. In the worst case it makes some tens of percents. Currently, such an accuracy of T 
calculation is most likely unattainable, even if not only the optical model of a nucleus, but 
also more modern nuclear models are at use. As a bright example, it is possible to indicate 
the comparison of calculated and experimental values of strength functions of (d,p) or (d,t) 
type reactions. Now it is fulfilled for many nuclei in experiments with high resolution at 
Tandem Accelerator of the University and Technical University of Munich up to the ap- 
proximately half of neutron binding energy. In all cases, the significant information on the 
excited levels structure, which allowed one to make a comparison with the nuclear theory 
with the excitation energy of up to several MeV, is obtained through these experiments. 
Such techniques as [5] do not take this fact into account. 

The comparison of (see, for example, [6]) the results obtained in Munich with calcula- 
tions within modern nuclear models shows that details of a fragmentation process of any 
states of nuclear potential over real levels with excitation energy up to several MeV can 
be reproduced only with a significant error. First of all, this circumstance is caused by an 
insufficient accuracy of the notions incorporated in the modern model on the nature of the 
nucleons interaction in a nucleus, as well as by the inevitable approximations of its account 
in anyone of them. The uncertainty of parameterization of the concrete nuclear models (op- 
tical, quasiparticle-phonon, models of interacting bosons and fermions, etc) always brings 
an additional error to the process of determining [5] p with the help of calculated T values. 
First of all, it is pointed at by the unavoidable and rather significant discrepancy between 
the experimental and model-calculated values of energies even for the most low-lying levels 
of the simplest structure. Moreover, it is also emphasized by the essentially greater one at 
the calculation of the nucleus parameters in case when excitation energy is higher than 1 - 
3 MeV. 

Discrepancy between the notions of the level density extraction from the nuclear evap- 
oration spectra and the technique of determining the excited levels structure by the levels 
excitation probability in nuclear reactions is a principal one. Basic thesis of a technique [6] 
follows from the experiment and cannot cause doubts. Reasoning from this, the existence of 
big systematical errors in T calculation is expected in case of using the theoretical models 
for the experimental determination of the level density [5] from the evaporation spectra. 

One should not exclude a possibility of another explanation of a discrepancy of p values 
determined by the first and the third techniques. For example, it may occur if a discrepancy 
of model notions and the experiment for spin dependence of the level density is not taken 



3 



into account at energies lower than several MeV. Particularly, as there is a possible, but the- 
oretically unconsidered, strong dependence of collective enhancement factor in level density 
[3] on a level's spin. If such an appreciable effect exists, it is displayed for spins J exceeding 
~ 6. It follows from maximal spin values of neutron resonances in 150 Sm, 177 Lu and from p 
value obtained for these and other nuclei. 

2. Very serious and, probably, unsolvable problems arise at determination [7,8] of p and k 
parameters of the cascade gamma-decay process of the nucleus from the gamma-rays spectra 
S of any nuclear reactions, in which the nucleus-product is excited up to the energy of 5-8 
MeV and higher. 

Rather simple modeling [9] of process of transferring errors 5S of the determination of the 
total gamma-spectra intensity into the ones of primary gamma-transitions [7], subsequently 
extracting the values of the level density and radiative strength functions from them, has 
shown that the value of 5p/p < 40 — 50% may be obtained within the framework of [8] type 
procedure. But it is possible only under the condition of all relative systematical errors of 
the measured gamma-spectra belonging to the region SS/S < 0.001 — 0.003 for all energies 
of emitted gamma-transitions and depopulated levels. This conclusion is just, at least, for 
the case of monotonous increasing or decreasing of the systematical errors SS as the nucleus 
excitation energy changes. Such an accuracy of measuring the intensity of total gamma- 
spectra by means of scintillation detectors is most likely unattainable for the present. It 
appears so, both because of the instability of their work, and of the impossibility of such an 
exact subtraction of Compton background from the instrumental spectra. Besides, in case 
when there is a nucleus with the high level density, the subtraction errors are the greatest. 

Additional and very essential systematical errors in the data for p and k, obtained by 
means of a technique of the primary gamma-transitions spectra determination [7] in corre- 
sponding experiments, are conditioned by its authors applying their own [8] algorithm for 
searching the maximum of likelihood function. It is used instead of the updates of Gauss 
method developed for completely similar cases and well tested by mathematicians. It is 
impossible to estimate the systematical errors related to it, as the method [8] can provide no 
information on position and value of its false maxima unlike the modern variants of Gauss 
method. 

The assumption stated in [8] deals with the equality of radiative strength functions of 
gamma-transitions of equal multipolarity and energy for depopulating levels with different 
energies. It contradicts the basic ideas of modern nuclear models (see, for example, [10]) and 
results in the unknown systematical errors occurring both in p and in k. 

Therefore, there is a necessity to develop the essentially new techniques for p and k de- 
termination. These should provide their extremely possible reliability. As mentioned above, 
such techniques cannot be created on the basis of any of the model notions for obtaining 
the experimental values of level density and radiative strength functions used both sepa- 
rately, and combined. That is why any of the newly developed techniques should be both 
model-independent, and determine the level density and radiative strength functions of the 
reaction's emission products. It should be simultaneously accompanied with providing the 
data accuracy of the modern experiment for the whole of the excitation energies range of the 
studied nucleus from the most suitable, for this purpose, reaction products spectra. This 
operation should be realized within the framework of the insufficient opportunities (for the 
guaranteed problem solution) and of essential limitations of a mathematical statistics algo- 
rithm developed for such purposes. Exactly at this stage of determination of the investigated 
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phenomenon parameters from strong correlating experimental data there appear the most 
serious systematical errors for both level density and radiative strength functions. It should 
be noted that in the mathematical expressions, describing the experimental spectra, the 
nonlinearity of the sought parameters connection is exactly responsible for the significant 
reduction of the area of pairs' parameters values with a hundred-percent correlation con- 
cerning a case of the linear equation systems. On this account, the degenerate systems of 
the nonlinear equations can have infinite number of solutions, all of which are in the limited 
intervals of their values. 



3 The potential of modern experiment to precisely de- 
termine the main parameters of a nucleus 

On one hand, proceeding from the base principles of mathematical statistics, the extraction 
of a number of unknown parameters from the data of an experiment demands a superfluous 
number of values to be determined in its course. On the other hand, the system of the 
sought parameters functional connections with the measured values has to provide both 
unambiguity of determination and minimum of errors of the sought values. Most likely, these 
two conditions cannot be satisfied in cases of obtaining the level density p from evaporation 
spectra [5] of any nuclear reactions or the radiative strength functions k both from total 
spectra [7] of gamma-rays following neutron radiative capture and from nuclear reactions at 
the charged particles beams [8]. 

Therefore, the experimental data of other type are necessary in order to get main nuclear 
parameters with higher accuracy than achieved at the present. The analysis [11] of the 
experimental data on cascade gamma-decay of one (or several) of the highly excited levels 
(compound states), fixed at energy B n , satisfies the listed above requirements to the maximal 
degree. First of all, these are the two-step cascade intensities 

Irrm = E E Y^jr = E , m nxi < y ~f (1) 

following the thermal neutron radiative capture, connecting a compound state A and a group 
of low-lying levels /, in function [12] of energies E 1 of their primary gamma-transition. Nu- 
cleus excitation energy (energy of the intermediate level of the cascade Ej) is unambiguously 
determined by energy E\. E^ — B n — E\. Cases of full absorption of the two-step cascade 
gamma-transitions energy, necessary for the experimental definition of 7 77 , are concentrated 
in the narrow peaks of the sum coincidence spectrum (Fig. 1) obtained with the help of 
ordinary HPGe-detectors. In the same experiment an inevitable background is measured 
with the maximal possible accuracy. 

Functional (1) depends both on ratio of partial and total radiative widths T of primary E\ 
and secondary Ei gamma-transitions of cascades between levels A, % and / and on number 
of levels n(m) = p x AE excited in different energy intervals AE. A degree of detailed 
elaboration of finding the form of energy dependence the level density and radiative strength 
functions is determined by the optimum interval's width of an averaging-out of cascades 
intensities. Technically, its limiting value equals the HPGe-detectors resolution. However, 
combined with the inevitable partial width fluctuations, real capabilities of the modern 
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detector and computer facilities limit AE interval width, over which nucleus excitation 
energy of 50 keV value (or a little larger) is distributed. 

Cascade transition type (dipole electric or magnetic) and an excited intermediate level's i 
spins and parity are unequivocally determined by known .P values of A and / levels. Practical 
absence of cascades between levels with | J\ — Jf\ > 2 excludes the necessity of accounting 
for transitions with higher multipolarities in analysis like [11]. For convenience of direct 
comparison of the obtained radiative strength functions with E1-, and Ml-transitions in 
nuclei with different mass A, it is appropriate to determine them in the following form: 

k = T M /(E% x A 2/3 x D x ) (2) 

and to use their ratios k(Ml)/k(El) for mutual normalization of data, experimentally mea- 
sured near B n . 

Intensity of cascades in expression (1) is proportional to a derivative dk/dE and, as a 
first approximation, it is inversely proportional to p. It provides the maximal sensitivity 
of the experiment in the range of p lowest values. So, in the field of excitation energy, 
the influence of the nuclear structure on parameters of studied nuclear reaction should be 
maximal. The essentially new connection type between the sought parameters, compared 
with the ordinary evaporation and gamma-spectra, provides smaller influence of correlation 
parameters on their real error. 

The specified experiment has two sources of the ordinary, but potentially rather signifi- 
cant, systematical errors. They are connected to a possibility: 

(a) Full compensation of a divergence between the experimental and calculated intensity 
values of the cascades with primary transitions E\ by a divergence of an intensity opposite 
sign of cascades with secondary transitions from the same energy intervals and 

(b) Presence of an excessive systematical error at the normalization of J 77 intensities. 
The first type error is related to the impossibility of instrument determination of the 

quanta ordering in the bulk of observed cascades. To full extent it is visible in analysis [13], 
using for extraction of the information on p and/or k the fact that the intensities directly 
observed in the experiment and calculated for different tested functional dependences of 
parameters are equal to their values for all cascade quanta energies. According to [9] this 
error leads to that the experimental spectra can be precisely reproduced by the level density 
and the radiation strength functions distinguished from the sought ones by some tens of 
times. The indicated error is practically reduced to zero only in the case, when the intensity 
of cascades in function of their primary transition energy is determined from the experimental 
spectra [12]. The corresponding algorithm uses the experimentally proved fact: even in nuclei 
with the maximal level density there are pairs of energetically resolved peaks, corresponding 
to the cascades with Ei < 3 — 5 MeV, which concentrate the main part of sum cascades 
intensity. A threshold of their intensity registration is much less than the average value even 
when rather ordinary detectors are used. Cascades with primary gamma-transitions from 
the same energy interval are registered mainly in the form of continuous distribution of small 
amplitude. Total background in any experimental spectra of two-step cascades intensities 
is equal to a determination error of a substrate under peaks of full capture of the cascade 
energy by the HPGe-detectors (Fig. 1). Practically, very exact algorithm, realized in [14] 
in combination with the results obtained by nuclear spectroscopy, allows one to determine 
the quanta ordering with high accuracy. Consequently, it makes possible for one to find the 
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form of energy dependences of experimental values I^Ei) with a relative systematical error 
less than [15] 0.01 - 0.05. 

The approximation of distribution of the resolved cascades intensity random values in 
small energy intervals of their intermediate levels allows [15] both determination of the cor- 
responding error 5/ 77 , and direct estimation of the most possible level density in any energy 
intervals of cascade intermediate levels below pa 0.5B n . It happens within the framework of a 
hypothesis about Porter-Thomas's law on fluctuations of the primary transitions intensities 
with the subsequent extrapolation of the approximated distribution below gamma-quantum 
registration threshold. In all of the analyzed nuclei the level density either practically co- 
incides with the data [11], or differs from the latter much less, than that measured with 
techniques [5] and [8]. 

Practically, expected systematical error 5/ 77 of the intensity cascades used in a technique 
[11] is almost fully determined by an error of measuring of the most intense primary gamma- 
transitions of cascades with low-lying intermediate levels. Comparison of the corresponding 
data from [16] and [17] shows that its value may be estimated as pa 10 — 20%. Coefficients of 
an error 5/ 77 transfer onto the values Sp and Sk are estimated in [18] for even-odd isotopes 
of tungsten and osmium. It was made for all of their excitation energy values. An error 
variation in an interval of —25% < 5/ 77 // 77 < 25% changes values, found for the level 
density below 0.5B n , less than by 2 times, and for radiative strength functions - less than 
by 3 times. 

It should be noted that the error in determining of J 77 can be larger than ±25% in the 
nuclei where the thermal neutron capture cross section was measured with the big relative 
uncertainty. In these cases, the amplitude of the expected relative errors of level density and 
radiative strength functions can significantly increase as compared with the results obtained 
in [18]. For some nuclei [6], this problem can be completely removed by comparison of the 
measured gamma-transition intensities with the known absolute intensities of gamma-rays 
related to the following /3-decay of the compound nucleus. One can only assume, however, 
that the relative errors of the measured capture cross sections for the most nuclei represented 
below are, in average, close to zero and have different sign. 

The third systematical error, by decreasing of the importance, is introduced into the 
results of technique [11] by using the notions of the identical form of energy dependences 
of radiation widths of the primary and secondary gamma-transitions of given multipolarity 
and energy. For the first time, the problem was formulated in the obvious terms and partly 
solved [19] at the experimental p and k determination from intensities J 77 of the two-step 
cascades. Its partial solution is achieved by using the experimental information on a full 
cascade population of the maximal number of levels in the bottom part of the investigated 
excitation energies region. The main result of the joint analysis [20] states that the level 
density decreases at the account of radiative strength function dependence on the energy 
of depopulating level as compared with the data [11]. It regards the two-step cascades 
intensities distribution onto low-lying levels and the data on population the high-lying ones 
for two tens of nuclei. 

In a technique [11] the examined systematical error displays itself to much smaller degree 
than in [8]. Due to that the revealed in radiative strength functions of secondary gamma- 
transitions changes have a different sign for various excitation energies, they are summarized 
in Ti with a multiplier (2) and, consequently, are relatively less in changing the total 
radiative width of cascades intermediate level. This conclusion is just only for the obtained 



7 



by the present time data on intensities of two-step cascades to final levels with energy 
E f < 0.5 - 1 MeV. 

As in any other experiment, at the extraction [11,20] of the cascade gamma-decay pa- 
rameters, the existing ideas are used as well as constants and earlier established concrete 
nucleus's parameters. First of all, it concerns: 

- the number of low-lying levels determined by nuclear spectroscopy [21], 

- an average spacing D\ between the neutron resonances, T\ total radiative width [22] of 
a decaying compound state. 

The corresponding data bring also additional, but correlating for various techniques, 
errors into the found p and k values. They can be significant only if the existing ideas of a 
nucleus contain essential errors. For example, it is so when the probability of excitation by a 
neutron of a level with a given J n higher than neutron binding energy, strongly depends on 
its structure. In such a hypothetical case, value of D\ can be essentially overestimated. It 
may be also possible that the branching coefficients at the decay of even rather long-living 
levels noticeably depend on their excitation way. 

If not to take into account such exotic opportunities, the information on cascades of 
gamma-quanta allows one to determine p and k without using the model-calculated values 
with the smallest possible systematical error for any stable target nucleus at any neutron 
beam. It is also needed to carry out [23] an independent check of the found values of the 
level density and strength functions by comparing the calculated and experimental total 
gamma-spectra. It can be fulfilled at capture of both thermal and, in principle, resonance 
neutrons. 

Thus, comparison of possible systematical errors of p and k values in experiments of 
various types shows a basic advantage of techniques [11,20] over the known ways of deter- 
mining only the level density [5] or of simultaneous determining the same value together 
with absolute ones of radiation strength functions [8]. It is caused by the following: 

(a) due to the instrument selecting only cases of full capture of the cascade energy, 
high resolution and stability of semi-conductor detectors the two-step cascade spectra are 
measured with a smaller systematical error as compared with the gamma-spectra, used by 
technique [8]; 

(b) the value J 77 depends on absolute value of the level density (used by the alternative 
techniques in their existing [5,8] variants, spectra do not depend on the absolute p and k 
values) ; 

(c) the transfer coefficients of a spectrum's error onto the 5p and 5k values in method 
[11,20] is ~ 2 orders smaller than [8]; 

(d) owing to the condition (b) there is a physically determined and essential limitation of 
the interval of possible p and k values, precisely reproducing both the intensity of cascades in 
energy functions of their primary transition, and other functionals of cascade gamma-decay 
process of any compound nucleus. 

(e) especially, it is necessary to emphasize that fixing of spins levels, for which the level 
density is determined, is practically unequivocal in expression (1) and there were no similar 
experimental results on this parameter till now. 

It should be noted that there is an essentially unavoidable error of p and k extraction 
from cascade intensities. It is caused by an excess of the unknown parameters number in 
(1) over the number of the experimentally measured values J 77 . Accordingly, the concrete 
values J 77 can be reproduced with an equal x 2 by infinite set of parameters determining 
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them. But the region of possible p and k values is always limited, providing that the 
relation of strength functions of primary and secondary transitions of equal energy and 
multipolarity is unequivocally fixed for any energies of decaying levels. Its width is minimal 
in the case of using all of the available experimental information on the investigated nucleus 
for determining the cascade gamma-decay parameters. At following these conditions, an 
interval of the most probable p and k variation does not in reality exceed several tens of 
percents if the value J 77 , used in [11,20], is [12] function of energy of the only primary cascade 
transition. 

Potentially, obtained p and k values can be distorted because of the influence of all three 
cascade levels structure on their intensity. The structure of a decaying compound state 
can influence 7 77 in a wide excitation energy interval, intermediate - locally. The influence 
degree of the structure of two-step cascade final level on its intensity is evidently shown in 
[19,20]. Corresponding effects for cascades primary gamma-transitions can be estimated and 
reduced only at their being studied in many neutron resonances. At present, a variation of 
the level density in concrete nucleus with respect to its average general trend (Fig. 2-8) can 
be accepted as the top value of an influence effect of resonance structure on the determined 
p and k values. 

4 On a basic possibility of the precise model- 
description of the modern data on the experimental 
level density 

The general form of dependence of the most reliable modern values of the level density 
on the investigated nucleus excitation energy points at [11,20] the presence of at least two 
"step-like" structures below neutron binding energy with its faster increase in between than 
it is predicted by a notion of a nucleus, as a system of non-interacting Fermi-particles. It 
means that, at least, at two excitation energies the abrupt change in the excited levels wave 
functions structure occurs in a nucleus. The unique factor, known by the present to be 
capable of providing such a change is the breaking of nucleons Cooper pairs with addition 
of two quasiparticles to the existing ones, as well as fast increasing the level density at 
increasing excitation energy. 

In modern theoretical notions [3], the level density at a given nuclear excitation energy U, 
spin J and parity n is expressed through the density p qp only of the quasiparticle excitations 
and its vibration and rotational (for the deformed nuclei) enhancement coefficients K v n, r and 
K rot , respectively: 

p(U, J, tt) = p gp (U, J, ir)K vibr (U, J, 7r)K rot (U, J, n) = p qp (U, J, ix)K coU {U, J, vr). (3) 

For the further analysis of the experimental data it is expedient to unite coefficients of 
vibration and rotational increase in the level density in the general coefficient of its collective 
enhancement K co \\. In the level density from the analysis [11,20] the basic contribution to 
its value is brought by the effect of vibration. In an examined case the effect of rotational 
enhancement of the level density for the deformed compound even-odd nucleus is less than 
the experimental data error. Probably, it is a little more for nuclei with neutron resonance 
spins of J > 2. By the order of magnitude in a neutron binding energy range, it is expected 
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that K co u value for the complete level density is [3] in an interval: 10 < K co u < 100. There 
is no experimental information on dependence of K vibr on U, J, and n. Modern theoretical 
ideas of this account admit a significant change in K vibr when changes U up to the change 
of [24] its functional forms dependences on nucleus excitation energy. The presence [11,20] 
of rather reliable experimental data for the sums of radiation strength functions of dipole 
cascade transitions allows, basically, to solve this problem by creation of precise models of 
strength functions. This opportunity is caused (see, for example, [10]) by known distinctions 
in values of partial radiation width from a ratio of quasiparticle and vibration components 
in the structure of excited (decaying) levels. However, currently, there are no theoretical 
models of such level [3]. 

Therefore, the further analysis of the level density is possible to be carried out only within 
the framework of zero assumptions of K co u independence on the nucleus excitation energy 
in the interval of ~ 0.5 — 3 MeV up to B n . For the first time, it allows one to receive direct 
experimental information on partial density of quasiparticle excitations with various number 
of quasiparticles in the specified excitation energy intervals for some nuclei with their masses 
of 40 < A < 200. 

The possibility of determining the partial level density p n with a given number n of the 
excited quasiparticles for U nucleus energy excitation 



has been found by Strutinsky [25]. For the first time, he has obtained simple functional 
dependence of the nucleus excited states density (the second coefficient in (4)). In more 
general form than expression (4), the model takes into account the existence of proton and 
neutron Cooper pairs (in the first turn [2] - for fissionable nuclei). All kinds of its modern 
notions are included in reviews [3], being quoted from original publications. Modern state 
of methods of partial density calculation is analyzed by Betak and Hodgson in [26]. 

Practically, for comparison with the experiment [11,20] within the framework of existing 
[3] theoretical notions of a nucleus, it is necessary to choose value of the spin cutoff factor 
a = f(n, U) for the given Cooper pair and excitation energy, together with the energy U — E n 
of the excited quasiparticles. Density of single-particle levels g for the presented here nuclei 
is known from the data on neutron resonances. The parameter g can take into account 
shell inhomogeneities of single-particle spectrum and depend on U. These rather theoretical 
opportunities have not been taken into account in the analysis because of the inevitable 
increase in its conclusions uncertainty. 

4.1 Fitting conditions 

In first of the tested by us variants, the E n = 0.25g(A^ — A 2 ) + Se n , functional dependence 
suggested by A.V.Ignatjuk and Yu.V.Sokolov (see [24]) has been used. It is based on the 
idea of existence of 0.25(/Aq condensation energy in a nucleus at its transition from normal 
to superfluid state. The maximal number of decayed pairs in this variant is limited to value 
N = 5, basically, by possibilities of used fitting algorithm. The essential increase in N 
could reduce xj value at least up to equaling it to xl value of the second variant. But 
it may be necessary to insignificantly change the obtained conclusions. For the present 



(J + l)exp(-(j + l/2) 2 /(2a 2 )) g n (U - E n ) 



(4) 
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the parameter Se n for next decayed pair is determined only [2] in modeling. However, 
for comparison of (4) with the experiment it was selected so that to provide the best p 
fitting for the experimental values [20] (when they are absent - the data on level density 
from [11] are used). Correlation function A of the ground state of the even-even nuclei is 
accepted as equal to the experimental (determined from atomic masses) pairing energy of 
last neutron; whereas in odd-neutron nuclei A = 12.8/v^4 approximation was used. Energy 
dependence of correlation function A n for the Cooper pair of number N at nucleus excitation 
energy U was set within the framework of its theoretically determined and presented in [2] 
approximation. It assumes that n-quasiparticle (n = 2N) excitations can exist above the 
threshold energy U th . It may be defined by the following expressions: 

3.144(n/ra c ) - 1.234(n/n c ) 2 for n/n c < 0.424 (5) 

1 + 0.617(n/n c ) 2 for n/n c > 0.424 (6) 

Here, C = 0.25gA^ is the condensation energy, and n c = 0.791gA is the number of 
excited quasi-particles in the vicinity of the phase-transition point from the superfluid to 
normal state. 

The energy dependence of the pairing-gap parameter A n , that is required for calculations 
of the level density (4), can be parameterized in the form: 

^ = 0.996 - 2.3Q(n/n c ) 1 - 57 /(U/C) - 76 for U/C > 1.03 + 2.07(n/n c ) 2 ' 91 (7) 
Ao 

— — = otherwise (8) 

A 

To correctly account the spin dependence of the level density is a serious problem because 
of the ambiguity of [3] model notions of the spin cutoff factor. In both described here variants, 
the functional dependence [27] suggested by Fu, has been used. The corresponding function, 
written in language FORTRAN77, is taken from file RIPL-2 [3] with necessary parameters. 
As the breaking threshold for number N pair was selected only by comparing values of 
expression (4) with the experimental density of levels, so its value, suggested in [27], has 
been replaced with the best one obtained by us. Values used in calculation of constants 
g and the best fitting parameters En (for En = U t h) of expression (4), found in the first 
variant, are presented in tables 1 and 2. 

First of all, it is found out that in this representation the experimental level density 
requires accounting for not less than 5 partial densities for its reproduction. According to 
the notions [25], the effect of pairing takes a share of nucleus excitation energy, which is equal 
to 2Ao- For five breaking nucleon pairs the total energy of pairing is equal approximately to 
10 MeV for the heaviest of the nuclei that are included in the analysis. This energy is more 
than the neutron binding energy and, consequently, notions [24,28] together with the data 
such as obtained in [5] seriously disagree with our experimental data. 

The substantial problem is also presented by discrepancy of energy thresholds of various 
pairs' breaking. As a rule, for the investigated nuclei this threshold for the 5-th pair is 
significantly less than for the 4-th. Sometimes inversion is observed for other pairs as well. 
In some of the nuclei still more inversion of disintegrated pairs' thresholds is observed. The 
trivial explanation can consist of a divergence of model notions [2] and the experiment in 



Uth 
C 
Uth 

c 
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close vicinity to the quasiparticles pair-production thresholds. If it is so, results of the 
executed approximation of the experimental level density in the first variant analysis show 
that with the error of about 1 MeV of breaking energy of 3, 4 and 5-th nucleon pairs 
practically coincide. Within the framework of notions [2] about the form of Cooper pairs 
correlation functions energy dependence of excited nucleus, a picture, corresponding to the 
basic statement [28], is observable. This statement deals with the generalized model of a 
superfluid nucleus - its phase-transition between superfluid and normal states (but at the 
essentially smaller energy of such transitions). 

Model approximation of correlation function A n , presented in [2], is obtained on the base 
of the experimental data such as [5]. Their basic distinctive feature is that the speed of 
increase in the level density is smaller as compared to data of [11,20]. It is true at least for 
excitation energy some higher than 0.5B n . In this variant of model notions, at the increase 
in excitation energy the smooth enough change in A n can provide the greater dp/dU value 
only when the expression (4) accounts for five and more breaking pairs. An alternative 
opportunity consists in using other ideas of nucleon correlation functions of Cooper pairs in 
the excited nucleus. 

The smaller values of n can be obtained only with the more rapid decrease in A n , than it is 
predicts by eq. (7), as increasing of the excitation energy U. And equation (4) cannot present 
any other opportunity. The number of variants of the functional dependences satisfying this 
condition is great. In the second variant of the analysis for E n the following functional 
dependence has been used: 

U-E n = U- A ln[(U - U th )/(pA )}. (9) 

The only ground for using function (9) for this purpose is a logarithmic dependence of 
macrosystem thermal capacity in a point of second-order phase transition on its tempera- 
ture. But such dependence is shown only in an ideal case. In case of a mixture of helium 
isotopes, for example, the maximal thermal capacity decreases, whereas the degree of change 
increases at the increase of 3 He concentration. Therefore, dependence (9) can be accepted 
as the utmost possible A n estimation for pair number N at the energy U with an additional 
condition [25] that maximum E n value is in E n < nA /2. Probably, following this condition 
for the majority of nuclei gives overestimated E n values. Only for 74 Ge, 185 W, 192 Ir and 
196 Pt, it is required to increase the maximal value of E n in (9) by 1.1 - 1.5 times in order to 
obtain the minimum possible xl- Such correction is essential only for the second breaking 
Cooper pair. If presence of systematical errors in p is taken into account, the specified E n 
increase does not indicate its excessive divergence in various nuclei. It is also true for the 
known significant fluctuations of pairing energy of last nucleons pair. 

The best value of parameter p for all tested nuclei is about 2.2-2.3. Therefore, in the 
second of the investigated variants, the pairing influences the level density p qp (U, J,ir) for 
next breaking Cooper pair only in essentially limited, as compared to (7), energy interval 
Uth < U < U t h +pA . I. e., it is practically equal to the known value of gap width in the 
low-lying states spectrum of an even-even nucleus. 

In the carried out variants of the analysis, the coefficient K co u has been accepted as 
independent on excitation energy U, spin and parity of levels and unchanging at the increase 
of nucleus excitation energy. Its absolute value is almost completely determined by a ratio 
of the experimental level density and the density of two or three quasiparticle excitations. 
In a case, when this value for any number N > 1 Cooper pairs differs by some times, such 
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discrepancy is easily compensated by changes in U th within the limits of, several hundreds 
keV maximally. 

At the analysis of the experimental data in even-odd compound nucleus some ambiguity 
arises at a choice of threshold of 3-quasiparticle level excitation energy. If the threshold of 
their production is accepted as sufficiently high (about several MeV), in such nuclei, the 
required collective enhancement factor of one-quasiparticle level density will exceed similar 
value for even-even and odd-odd nuclei by a factor of some tens. Therefore, it has been 
postulated, that in the excitation energy interval of about 2-3 MeV above the ground state 
of even-odd nucleus, the base density is the one of three-quasiparticle excitations instead of 
the density of one-quasiparticle ones. 

In odd-odd compound nuclei the density of 2-quasiparticle excitations assumes that there 
is an excitation both of neutron, and proton quasiparticles. The type of quasiparticles for 
4 and more quasiparticle excitations in these nuclei, as well as in even-even ones cannot be 
established by the analysis carried out here. 

4.2 Results and their discussion 

Examples of the best approximation of the experimental data [11,20] for 4 nuclei with various 
parity of neutrons and protons are displayed in Fig. 2. It also shows the partial level density 
obtained in the first variant of the analysis. In Figs. 3-8 the similar data are presented 
for the most of the analyzed nuclei with partial densities obtained in the second variant of 
analysis. In tab. 2 the parameters of the approximated partial level density for the second 
variant of accounting for nucleon pairing interactions in nucleus are presented. Thresholds 
of excitation energy U t h, necessary for calculation by expressions (4) and (9), are averaged 
in it after division of En onto the approximated value of correlation function Ao = 12.8/v^4 
for nuclei with different neutron and proton parity, separately. 
From these data quite unequivocally follows, that: 

1. The first step-like structure in level density [11,20] is caused by existence of, at least, 
two quasiparticles in nucleus of any type; 

2. For the precise reproduction of the level density (comparable with an experimental 
data error) it is required to postulate the breaking from three up to five and more such pairs; 

3. 2-3 and more nucleon Cooper pairs can have practically equal breaking threshold if only 
the effect of nucleons pairing of any pair demonstrates itself in a wide [2] interval of nucleus 
excitation energy. If not to take into account the circumstance that the approximation of 
experimental data on p both from [11] and [20] necessarily inverts thresholds of breaking at 
least of the fourth and fifth pairs, the proximity of their values allows one to speak about 
practically observable simultaneous breaking of several Cooper pairs. That is, about [28] a 
nucleus phase-transition from superfluid to the normal state; 

4. The coefficient K co u, taken from the data [20], has practically equal value in both 
even-even and even-odd compound nuclei and is significantly higher in odd-odd ones (for 
the level density from a variant of the analysis [11]). This difference is qualitatively explained 
by value and error distinction sign of the level density given by [11] as regard to [20]; 

5. From the data of approximation, presented on Figs. 3-8, it is visible, that in the 
second variant of the analysis the experimental level density can be well reproduced at 
the account of only 6 or 7-quasiparticle excitations. Moreover, representation of correlation 
functions as (9) allows one to obtain small enough value of the level density of n-quasiparticle 



13 



excitations, having lower energy than U t h- So, the known from nuclear spectroscopy fact of 
the 2-quasiparticles excited levels presence within the limits of nucleus excitation energy 
somewhat lower than 2A is thus explained; 

6. In Tables 1 and 2 data on comparison of average values of collective enhancement factor 
of the level density are collected and show smaller scatter in the second variant of the analysis. 
In this variant differences < E 2 /A > - < E 1 /A >^ 4 and < E 3 /A > - < E 2 /A >k, 2 
coincide in nuclei of various types within the limits of an error; 

7. In [25] it is predicted that the level density of odd A nucleus corresponds to that of 
even nucleus with excitation energy U + A , accordingly. Results of the first variant of the 
analysis (Tab. 1) do not prove this prediction to be true. However, in the second one it is 
carried out within the limits of E n energies determination error; 

8. The shell effects demonstrate itself to the maximum extent in near magic nuclei being 
close to iV=82 and 126. Their exhaustive reproduction is impossible within the framework 
of expression (4). More exact approximation can demand to revise the basic assumption [25] 
about equidistant character of single-particle spectrum of near magic nucleus and to account 
[24] for shell inhomogeneities of single-particle spectrum. It may also be required to account 
for change in K coU value for different nuclei and excitation energies, not taken into account 
by performed analysis. 

Therefore, it is undoubtedly required to obtain better experimental data on gamma- 
gamma coincidences and further development of algorithms of their analysis. First of all it 
concerns the determination of E n energy values, connected with collective type excitations 
near a breaking threshold of pair number N. The problem of K co u value determination at a 
given excitation energy U is directly connected to this circumstance. For even-even nuclei it 
is visible from the given tables and figures. Approximation of p is done in them starting with 
U ~ 2 MeV. Just this is the explanation of the discrepancy of E 1 energy value with 2A . 
As long as in even-even nuclei of A region there are only levels of vibrational type, so far it 
is necessary to accept K co u as oo in this excitation region. The available experimental data, 
unfortunately, do not give an opportunity to estimate this problem near to values E 2 and E 3 
basically because of the lack of theoretical representations on this account and mathematical 
methods of the analysis of the experimental data. 

Despite of the specified ambiguity, volume and quality of the obtained [11,20] experimen- 
tal data presents extremely favorable possibilities for the greatest possible development of 
modern models of the level density. The obtained values of the most probable parameters 
of the experimental data approximations can be used for p calculation in nuclei that are 
relatively far from the magic ones. It may be done either directly or at interpolation of the 
data given by Tables 1 and 2, or by using the corresponding average values. It is certainly 
necessary to develop new model ideas of the nucleon correlation functions energy dependence 
for near magic nuclei in the whole region of E ex < B n and about p vibration enhancement 
factor. 

It is also necessary to reveal the degree with which the initial levels A structure influences 
on the two-step cascade intensity and to estimate the p and k systematical error connected 
to this circumstance. 
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5 Conclusion 



There is an urgent need for both the experimental and the theoretical determination of the 
excited level density in a given nuclear reaction, as well as the probability of its products 
emission for all possible excitation energies with a guaranteed systematical error not ex- 
ceeding several tens of percents. That is, it is to be accomplished with the same accuracy 
of measuring various spectra of nuclear reactions and their cross sections, as carried out, 
presently. Only such reduction of available systematical errors of the experimental p and k 
determination allows one to observe dynamics of nucleus structure change at its excitation 
energy variation and to reveal factors causing it. There is no other possibility to develop new 
nucleus model and to specify the existing ones. Known techniques [5,7,8] of p and k determi- 
nation are incapable of solving this problem. It is true either because of the impossibility of 
getting accuracy-specified model-designed transmission coefficient for nuclear reaction's nu- 
cleon products, or because of practical inaccessibility of the required measurement accuracy 
of gamma-rays spectra, accompanying nuclear reactions with the charged particles emission. 

Specific feature of the model-free simultaneous extraction [8,11,20] of both the level den- 
sity and strength functions of emission of the nuclear reaction's registered products is nec- 
essary to determine the intervals of their values, reproducing experimental spectra with the 
required accuracy. Essentially, it is provided for by the irremovable ambiguity of such prob- 
lem's solution both for [8] and [11,20] techniques. The marked circumstance is connected 
with degeneracy of the equation systems, connecting values of the sought parameters with 
the measured spectra intensities. Model notions and hypotheses (first of all, on independence 
of strength functions on a nucleus's excited levels structure) reduce or fully eliminate [5] the 
equation systems degeneracy. But they necessarily lead to an unknown systematical error p 
and k. Therefore, reliable results for p and k can be only obtained by using a mathemati- 
cally correct technique of the experimental data processing. The data should most precisely 
reflect this circumstance. 

The accumulated, up to the present, data file on two-step cascades following thermal 
neutron radiative capture, points to the existence of systematical discrepancy in experimen- 
tally obtained data for p and k and their model-calculated values. These data also specify 
excessive idealization of notions about a nucleus, used by techniques of [5,7,8] type for the 
extraction of its major parameters. Presence of significant systematical discrepancy between 
the experiment and theory allows one to hope on the possibility of a substantial improvement 
in model notions about properties of a nucleus at its excitation being lower than nucleon 
binding energy. 
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Table 1. 

The best values of the collective enhancement factor K coU for analyzed nu- 
clei, En is the value of breaking thresholds of pair number N, MeV. 
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Table 2. 

A J is the interval of spins for which the experimental value p is determined, parameter p of 
expressions (9) and the best values of the collective enhancement factor K co u for analyzed 
nuclei in the second variant of the analysis. E N , MeV is the value of breaking thresholds of 
pair number N. R — xt/Xf ~ relations of criteria of fitting quality of the second (s) and the 
first (f) variants of the analysis. 
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Fig. 1. It presents the main part of the sum coincidence spectrum for the target enriched 
in 124 Te. Full energy peaks are labeled with the energy (in keV) of final cascade levels. 
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Fig. 2. It shows examples of the experimental data on approximation for nuclei 125 Te, 128 1, 
168 Er and 177 Lu by the partial level density in the first variant of the analysis. Full points 
with error bars represent the experimental data [20], open points show the data [11]. A thin 
dotted line indicates partial density points display their sum. A solid line presents the sum 
of the partial level density from the second variant of the analysis. Thin line stands for the 
level density calculated within model [29]. 
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Fig. 3. The same as in Fig. 2 for nuclei: 40 K, 60 Co, 71 ' 74 Ge, 80 Br, 114 Cd. Here, the partial 
level densities are shown as obtained from the second variant of the analysis. 
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Fig. 6. The same as in Fig. 2 for nuclei: 166 Ho, 170 Tm, 174 Yb, 176 Lu, 181 Hf, 182 Ta. 
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Fig. 7. The same as in Fig. 2 for nuclei: i83,i84,i85,i87 W; i88,i90 Os . 
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